NASA TM X-1339 


COMPUTER PROGRAM FOR THERMODYNAMIC PERFORMANCE 

OF BRA YTON- CYCLE SPACE -POWER SYSTEMS 

By Arthur J. Glassman 

Lewis Research Center 
Cleveland, Ohio 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


For sole by the Clearinghouse for Federal Scientific and Technical Informoticn 
Springfield, Virginia 22151 ■ 'rice $1.00 








; 













J .• ,» t 4* 

' . -v, •'■/■'■ 





COMPUTER PROGRAM FOR THERMODYNAMIC PERFORMANCE 
OF BRAYTON-CYCLE SPACE-POWER SYSTEMS 
by Arthur J. Glassman 
Lewis Research Center 

SUMMARY 

A computer program used for the calculation of the thermodynamic performance of 
one- and two-shaft Brayton-cycle space-power systems is presented. Systems that can 
be analyzed include those with and without reheating and/or intercooling. Provision is 
made in the program to account for turbine-coolant flow. 

Inputs required for the program include die component performance parameters and 
the cycle temperature variables. Performance outputs from the program include cycle 
efficiency and prime radiator area. Sample input and output are presented to demonstrate 
the use of the program. 


INTRODUCTION 

One of the first phases in the analysis of any power cycle is the determination of the 
thermodynamic performance parameters and optimum cycle- variable values pertinent to 
the particular application. Such analyses have been performed for Brayton-cycle space- 
power systems to determine the effects of cycle temperature variables, component per- 
formance parameters, reheating and/or intercooling, and turbine-coolant bypass flow on 
cycle efficiency and prime radiator area. Examples of these analyses are found in refer- 
ences I and 2. These cycle studies are performed most rapidly and conveniently with a 
computer. 

The computer program used to obtain the thermodynamic performance results in 
references 1 and $and associated studies is presented herein. The analysis equations 
comprising the program are given in these references. One- and two- shaft systems with 
and without reheating and/or intercooling can he analyzed. It is also possible to account 
for die effects of turbine-coolant flow. The features of this program are discussed and 



the FORTRAN listing is presented. Sample input and output are included to demonstrate 
use of the program. 


PROGRAM DESCRIPTION 

The computer program calculates cycle thermodynamic performance for one- and 
two-shaft Brayton-cycle space-power systems with and without reheating and/or inter- 
cooling. A schematic diagram of the general system being analyzed along with calculation 
station identification numbers is presented in figure 1. The system variations that can be 
analyzed are obtained by eliminating components such as the reheater, intercooler, and 
high-pressure turbomachinery through input specifications. The thermodynamic analysis 
equations that comprise the program were derived in references 1 and 2. The general 
features of the program are discussed herein and detailed descriptions of the input and 
output are presented. The FORTRAN program and associated variables are listed in the 
appendixes. 


General Features 


Basic systems . - The program is set up to analyze five basic systems, four of which 
are included in the general system shown in figure 1. These four systems, in which the 



Figure 1. - Brayton-cycle system with Intercooling and reheating. 
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alternator is on the low-pressure shaft, include those with (1) neither reheating nor inter- 
cooling, (2) reheating only, (3) intercooling only, and (4) both reheating and intercooling. 
In the fifth system, which also has both reheating and intercooling, the alternator is on 
the high-pressure shaft. As seen in figure 1, only one stage of reheating or intercooling 
is considered. The system to be analyzed is specified by the input variable KASE, as will 
be described in the section Input, and the program automatically selects the proper anal- 
ysis path. 

Shaft arrangements . - As seen in figure 1, the system is considered to have a two- 
shaft turbomachinery arrangement when reheating and/or intercooling are used. With 
neither reheating nor intercooling, either one-shaft or two-shaft arrangements can be 
analyzed. The selected shaft arrangement and particular turbine-work split for a two- 
shaft arrangement are specified by the input variable ST, which expresses the ratio of 
high-pressure-turbine work to low-pressure-turbine work. Appropriate selection of ST, 
as will be discussed in the section Input, yields shaft arrangements including a one-shaft 
system (the low-pressure shaft in figure 1), various two-shaft two-compressor systems 
(representing variations in shaft-work split), and a two-shaft one-compressor system 
(compressor on the high-pressure shaft). 

Turbine-coolant flow . - The effects of turbine-coolant flow are analyzed according to 
the bypass model specified in reference 2; that is, coolant flow is assumed to bypass the 
turbine and then mix with the main stream. Coolant flow, as seen from figure 1, can 
originate either from the compressor outlet (station 5) or the recuperator outlet (station 6) 
as desired. This coolant can then rejoin the main stream downstream of either turbine. 
The amount, source, and destination of the bypass flow are specified by the input vari- 
ables YIP, Y1PFR6, Y2, and Y2FR6 as described in the section Input. 

T emperature ratio variation. - For each set of input parameters, the two cycle tem- 
perature variables can be varied parametrically between specified limits. These vari- 
ables are cycle temperature ratio (T41), which is the ratio of compressor-inlet to 
turbine-inlet temperature, and turbine temperature ratio (T21), which is the ratio of 
turbine-exit to turbine-inlet temperature. The limits and increments for these tempera- 
ture ratios are specified by the input variables T41MIN, T41DEL, T41MAX, T21MIN, 
T21DEL, and T21MAX. For each value of cycle temperature ratio, the cycle computa- 
tions are made for the range of turbine temperature ratios, and then the maximum cycle 
efficiency within that range and the associated turbine temperature ratio are determined. 

Reheat and Intercool temperatures. - For any given case, the temperatures leaving 
the intercooler or reheater can be varied. These temperatures are expressed as ratios 
of the first-compressor-inlet or first-turbine-inlet temperatures, and are specified by 
the Input variables C4 and Cl. 

Program organization. - There is a main program, ETACY2, and three subroutines, 
RADTRA, MAXIM, and ZERO. Main program ETACY2 controls all the input and output 
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TABLE L - INPUT FORM WITH SAMPLE DATA 

























and computes all temperatures and turbomachinery pressure ratios, as well as cycle ef- 
ficiency and weight flow, for the cycle. Subroutine RADTRA computes the prime radiator 
areas required for both primary heat rejection and intercooling. Subroutine MAXIM pro- 
vides the logic for determining the maximum cycle efficiency and associated turbine tem- 
perature ratio for each cycle temperature ratio. Subroutine ZERO provides the logic for 
determining the turbine work split that yields a two-shaft one-compressor arrangement. 
These programs are listed and the variables identified in the appendixes. 


Input 

The input values for this program include component efficiencies, cycle loss pressure 
ratio, gas specific heat ratio, turbine work split, turbine-inlet temperature, space-sink 
temperature, radiator emissivity, gas heat-transfer coefficient in the radiator, 
reheat temperature, intercool temperature, cycle temperature-ratio range and increment, 
turbine temperature-ratio range and increment, and integers specifying the appropriate 
system and giving instructions for further input. 

In table I, sample data are presented in the required form for input. A field width of 
10 columns is allowed for all variables except those on the first card, which contains only 
2 one-digit integers. If there is no decimal point included in the field, the four rightmost 
elements of each field are taken to represent four decimal places. The input variables 
are defined in the following list: 

KASE specifies system as follows: 

1 - intercooling and reheating, alternator on high-pressure shaft 

2 - intercooling only, alternator on low-pressure shaft 

3 - intercooling and reheating, alternator on low-pressure shaft 

4 - no intercooling or reheating, alternator on low-pressure shaft 

5 - reheating only, alternator on low-pressure shaft 

KREAD specifies FORTRAN statement number (1, 2, or 3) to which program returns 
for more input data 

ETAT1 efficiency of turbine on high-pressure shaft, fraction 
ETAT2 efficiency of turbine on low-pressure shaft, fraction 
ETAC1 efficiency of compressor on low-pressure shaft, fraction 
ETAC2 efficiency of compressor on high-pressure shaft, fraction 
RL cycle loss pressure ratio 

E recuperator effectiveness, fraction 
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GAMMA 

T41MIN 

T41DEL 

T41MAX 

T21MIN 

T21DEL 

T21MAX 

Cl 

C4 


T1 

TS 

EPS 

HR 

DLMIN 

YIP 

Y1PFR6 

Y2 

Y2FR6 

DLT1PM 


gas specific heat ratio 

minimum value of cycle temperature ratio, (ratio of compressor-inlet to 
turbine-inlet temperature) 

cycle temperature- ratio increment 

maximum value of cycle temperature ratio 

minimum value of turbine temperature ratio (ratio of turbine-exit to turbine- 
inlet temperature) 

turbine temperature-ratio increment 

maximum value of turbine temperature ratio 

ratio of reheat temperature to turbine-inlet temperature 

ratio of intercool temperature to compressor-inlet temperature 

ratio of high-pressure-turbine work to low-pressure-turbine work (ST = 0 
yields a one-shaft arrangement and a negative ST yields a two-shaft, one- 
compressor arrangement) 

turbine-inlet temperature, °R 

space-sink temperature, °R 

radiator surface emissivity 

radiator gas heat-transfer coefficient, Btu/(hr)(ft^ prime area)(°R) 

convergence increment for radiator area computation 

ratio of high-pressure-turbine coolant flow to compressor flow 

fraction of YIP originating from recuperator outlet 

ratio of low-pressure-turbine coolant flow to compressor flow 

fraction of Y2 originating from recuperator outlet 

convergence increment for temperature after mixing of main stream and 
coolant 


Output 

The printed output from the program Includes pertinent input values and the following 
computed values: cycle efficiency, prime radiator area, a weight flow parameter, turbo- 
m ac hiner y pressure ratios, and ratios of all temperatures around the cycle to turfaine- 
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OUTPUT FORM WITH SAMPLE DATA 
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inlet temperature. All temperatures and pressures are stagnation values. Prime radi- 
ator areas are computed on the basis that the gas cooler and intercooler are radiators. 

Sample output corresponding to the sample input is given in table n. The top line is 
a program identification title; the second line specifies the case being analyzed. Next are 
two sets of two lines each, consisting of pertinent input variables and their associated 
values. The printed names for all but one of these input variables, except for added paren- 
theses in some cases, corresponds exactly to those listed in the section Input The one 
exception is SPLIT (T) which corresponds to the input variable ST. The computed results 
are then presented in groups, each group being for a different cycle temperature ratio 
(T4/T1) starting from the lowest (T41MIN) and going to the highest (T41MAX). Three 
such groups are shown in table n for cycle temperature ratios of 0. 275, 0. 300, and 0. 325. 
The first line in each group identifies the cycle temperature ratio (T4/T1) and gives the 
ratios of intercool temperature to turbine-inlet temperature (T4PP/T1, which is a func- 
tion of the input variable C4) and reheat temperature to turbine inlet temperature 
(T1PP/T1, which corresponds to the input variable Cl). With no intercooling or reheat- 
ing, the pertinent temperature ratio is artificially set to zero for printout purposes only. 
Next is a line of headings which are identified in the following list. Then there is one line 
of results for ea»;i turbine temperature ratio (T2/T1). The last line of each group shows 
the values that correspond to maximum cycle efficiency for that cycle temperature ratio 
group. The column headings for each group are identified as follows: 


T2/T1 

ETACY 

WCT/P 


T1P/T1 

T3/T1 

T4P/T1 

T5/T1 

T6/T1 

PI/PIP 

P1PP/P2 

P4P/P4 

P5/P4PP 


ratio of turbine-exit to turbine-inlet temperature 
cycle efficiency, fraction 

weight flow parameter, equal to weight flow (lb/sec) times specific heat 
(Btu/(lb)(°R)) times turbine-inlet temperature (°R) divided by shaft power 
to alternator (kW) 

ratio of high-pressure-turbine-exit to turbine-inlet temperature 

ratio of gas cooler-inlet to turbine-inlet temperature 

ratio of low-pressure-compressor-exit to turbine-inlet temperature 

ratio of high-pressure-compressor-exit to turbine-inlet temperature 

ratio of gas-heater-inlet to turbine-inlet temperature 

pressure ratio across high-pressure turbine 

pressure ratio across low-pressure turbine 

pressure ratio across low-pressure compressor 

pressure ratio across high-pressure compressor 
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SPLIT (C) 

ARAD/P 

AINT/P 

ATOT/P 


ratio of low-pressure-compressor work to high-pressure-compressor work 
(artificially set to zero when ST = 0) 

o 

specific prime radiator area required for primary heat rejection, ft /kW 

2 

specific prime radiator area required for intercooling, ft /kW 

2 

total specific prime radiator area, ft /kW 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, November 15, 1966, 
128-31-02-25-22. 
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APPENDIX A 

MAIN PROGRAM ETACY2 

Main program ETACY2 controls all input and output, contains all analysis -path 
branching logic, and computes all temperature ratios around the cycle, turbomachinery 
pressure ratios, cycle efficiency, and weight flow parameter. The thermodynamic equa- 
tions used are those presented in appendixes B of references 1 and 2. 


PROGRAM VARIABLES 

All items marked input variable are identified in the section Input. The variables for 
ETACY2 are the following: 

AINOP specific prime radiator area for intercooling 

AROP specific prime radiator area for primary heat rejection 

ATOP total specific prime radiator area 

Cl input variable 

C4 input variable 

CFUNC function of COEFC 

2 

COEFB coefficient B of quadratic equation x 

COEFC coefficient C of quadratic equation 

DLMIN input variable 

DLT1PM input variable 

E input variable 

EE product of recuperator effectiveness and fraction of compressor flow that 

passes through recuperator 

EPS input variable 

ETACY cycle efficiency 

ETAC1 input variable 

ETAC2 input variable 

ETANUM numerator of expression used to compute ETACY 


+ Bx + C = 0 where x equals T4P1 
+ Bx + C = 0 where x equals T4P1 


10 


ETAT1 

ETAT2 

EXP 

GAMMA 

HR 

IN 

IND 

J 

KASE 

KREAD 

M 

P1PP2 

P1P1X 

PUP 

P21PPX 

P4P4 

P4P4X 

P54PP 

P54PPX 

RCX 

RL 

RTX 

SC 

ST 

STT 

STY 

TEST 


input variable 
input variable 
function of GAMMA 
input variable 
input variable 

logic indicator for determination of ST value that yields SC = 0 

logic indicator for determination of maximum value of ETACY 

logic indicator for determination of program path during ETACY maximization 
iteration 

input variable 

input variable 

logic indicator for determination of T1PM1 

pressure ratio across low-pressure turbine 

isentropic tempera tuTfe^atio across high-pressure turbine 

pressure ratio across high-pressure turbine 

isentropic temperature ratio across low-pressure turbine 

pressure ratio across low-pressure compressor 

isentropic temperature ratio across low-pressure compressor 

pressure ratio across high-pressure compressor 

isentropic temperature ratio across high-pressure compressor 

isentropic temperature ratio corresponding to total pressure ratio across both 
compressors 

input variable 

isentropic temperature ratio corresponding to total pressure ratio across 
both turbines 

ratio of low-pressure-compressor work to high-pressure-compressor work 
input variable 
input value of ST 

value of ST adjusted for coolant bypass, 
value tested for convergence of T1PM1 
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TS 

TS1 

TT1PM1 

T1 

T1PM1 

T1PP1 

TIPI 

T2M1 

T21 

T21BEG 

T21DEL 

T21END 

T21MAX 

T21MIN 

T31 

T4PP1 

T4P1 

T4P4 

T41 

T41DEL 

T41MAX 

T41MIN 

T51 

T54PP 

T61 

WPARAM 

XT1P1 

Y 


input variable 
ratio of TS to T1 
temporary storage for T1PM1 
input variable 

ratio of high-pressure-turbine-exit temperature, after mixing with coolant, 
to turbine- inlet temperature 

ratio of reheater- exit to turbine-inlet temperature 

ratio of high-pressure-turbine-exit to turbine-inlet temperature 

ratio of low-pressure-turbine- exit temperature, after mixing with coolant, 
to turbine - inlet- temperature 

ratio of low-pressure-turbine-exit to turbine-inlet temperature 
lowest value of T21 used in ETACY maximization iteration 
input variable 

highest value of T21 used in ETACY nuugjpization iteration 
input variable 
input variable 

ratio of gas cooler-inlet to turbine-inlet temperature 
ratio of intercooler- exit to turbine-inlet temperature 
ratio of intercooler-inlet to turbine-inlet temperature 
temperature ratio across low-pressure compressor 
ratio of compressor-inlet to turbine-inlet temperature 
input variable 
input variable 
input variable 

ratio of high-pressure-compressor-exit to turbine-inlet temperature 

temperature ratio across high-pressure compressor 

ratio of gas-heater-inlet to turbine-inlet temperature 

weight flow parameter (see section Output) 

storage for calculated TIPI s ^ ‘ / ' 

lower limit test value for TIPI 


12 



YTOT 


ratio of total coolant to compressor flow 
input variable 


YIP 

Y1PFR6 input variable 

Y2 input variable 

Y2FR6 input variable 

Y5 ratio of coolant flow from compressor exit to compressor flow 

Y51P ratio of high-pressure-turbine coolant flow from compressor exit to com- 
pressor flow 

Y52 ratio of low-pressure-turbine coolant flow from compressor exit to com- 

pressor flow 

Y6 ratio of coolant flow from recuperator exit to compressor flow 

Y61P ratio of high-pressure-turbine coolant flow from recuperator exit to com- 

pressor flow 

Y62 ratio of low-pressure-turbine coolant flow from recuperator exit to com- 

pressor flow 


PROGRAM LISTING 


BIBFTC ETACY2 

C BRAYTON CYCLE - TWO SPOOL - VARIABLE MORK SPLIT 
C CASE 1 - INTERCOOL AND REHEAT - ALT. ON HIGH PRESSURE SHAFT 

C CASE 2 - INTERCOOL ONLY - ALT. ON LOU PRESSURE SHAFT 

C CASE 3 - INTERCOOL AND REHEAT - ALT. ON LOU PRESSURE SHAFT 

C CASE 4 - NO INTERCOOL OR REHEAT - ALT. ON LOU PRESSURE SHAFT 

C CASE 5 - REHEAT ONLY - ALT. ON LOU PRESSURE SHAFT 

C RASE DENOTES THE SPECIFIED CASE 

C ENO OF PROGRAN TRANSFER TO STATEMENT NUMBER KREAD 

C IF INPUT ST IS NEGATIVE* SC HILL EQUAL ZERO 

C 

1 REAO (S* 9991 RASE » KREAD 

READ IS* 1000) ETAT1*ETAT2*ETAC1 *ETAC2 
REAO IS* 1000) RL*E*GANNA 
REAO 15*1000) T41MIN*T410EL*T41MAX 
REAO (5*1000) T21NIN*T210EL*T21MAX 
T 41NAX*T 4 1MAX+T41DEL* « 1 
T 21MAX«T2 1MAX+T210EL* • 1 

2 REAO 15*1000) C1*C4*ST 
IFtCl.EQ.O.O) GO TO 1 
STT-ST 

REAO (5*1000) T1*TS*EPS*HR*0LM1N 

3 REAO 15*10.00) V1P«Y1PFR6*V2*Y2FR6*0LT1PN 
IF! Y1P.EQ.1.0) GO TO 2 
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15 WRITE 16.1010) 

XT1P1*1.0 

GO TO (30* 20.40*41 . 42 ) »KASE 
20 WRITE (6,1021) 

GO TO 50 

30 WRITE (6.1022) 

GO TO 50 

40 WRITE (6.1023) 

GO TO 50 

41 WRITE (6.1024) 

GO TO 50 

42 WRITE (6,1025) 

50 WRITE (6,1030) ETAT1,ETAT2,ETAC1,ETAC2,RL,E, GAMMA, Cl, C4.STT 
WRITE (6.1035) T1,TS,EPS.HR,Y1P,YIPFR6,Y2,Y2FR6 
Y62*Y2FR6*Y2 

Y61P-Y1PFR6*YIP 

Y6=Y62*Y61P 

YT0T-Y2+YIP 

Y5*YT0T-V6 

Y52=Y2-Y62 

Y51P*YIP-Y61P 

EE*E*( 1.0-Y5I 

T41-T41MIN 

T21-T21MIN 

T51-T21 

T61-T21 

TTIPM1-0.0 

J*0 

M*0 

55 T21BEG-T21MIN 
T21EN0*T21MAX 
T4PPl»C4*T4l 
TIPPl-Cl 

IFIKASE.EQ.2.0R.KASE.EQ.4) T1PPI*0.0 
IF(KASE.EQ.4.0R.KASE.EQ.5I T4PP1-0.0 
WRITE (6.1045) T41.T4PP1.T1PP1 
49 IFISTT.LT.-.OOl) ST*1.0 
IN-l 

51 4F(KA$E«EQ.1.0R.KASE«EQ.3.0R.KA$E.EQ.5) GO TO 60 
STY«ST*( 1.0-Y2I/C l.O-YTOT) 

TlPMl-(( l.-YT0T)*(l.+STY*f21)*Y5lP*T51*Y6lP*T6l)/((l.-Y2)*(l.+$T{) 

IF( T1PM1.LT. T21) T1PM1-T21 

TlPl»l.-STY*(TlPMl-T21) 

T1PPI*TIPMI 

IF(YIP.EQ.O.O) TTIPM1-T1PM1 
TEST«A9S( TIPMI-TT1PH1) 

IFITEST.LE.DLTIPM) M*i 
TTIPMl-TIPMl 
GO TO 70 
60 T1PP1-C1 
N«1 

STV«ST*( l.O-Y 2) /( l.O-YTOT) 

T1P1*1.-STY4(TIPP1-T2U 

XT1P1*TIP1 

Y»l*-6TATl 

IF(TIP1«LE.Y) TIP l- Y+. 01 J 
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70 T4PPl»C4*T 41 

P1PIX»1„-1 i.-TIPlI/ETATi 
P21PPX»1.-1 1.-T21/T1PP l 1/ETAT2 
RTXsU/(PlPlX*P21PPX) 

EXP*J GAMMA-1.) /GAMMA 
RCX»RTX/Pl**EXP 
IF(KASE.EQ.l) GO TO 80 
IF1KASE.EQ.4.0R.KASE.EQ.5) GO TO 85 

75 T5l*< i.-YT0T)*ll.-TiPn+T4PPi 
P54PPX*1. ♦( T51/T4PP1- 1. )*E TAC2 
P4P4X*RCX/P54PPX 
IF1KASE.GE.4) GO TO 90 
T4P4*i.«-lP4P4X-l.)/ETACl 
T4Pl=T4P4*T4i 
GO TO 90 

80 T4P1*( 1.-Y2>*IT1PP1-T21>*T41 
P4P4X*1.+1 T4P 1/T41-1. )*ETAC1 
P54PPX*RCX/P4P4X 
T54PP*i«*l P54PPX- l. 1/ETAC2 
T51»T54PP*T4PP1 
GO TO 90 

85 COEFB* 1 1 .-ETAC l-RC X*E TAC 1*E TAC2* I 1 .-YTOT I * 1 1. -Tl PI ) /T4I I *T4 1/ET AC1 
COEFOT41M U-ETACl)*ETAC2*(l.-YT0T)*f 1.-T1PU/ETAC1 
CFUNO-SQRTt 4.*C0EFC I 
IF1C0EFB.GT.CFUNC) GO TO 400 
T 4Pi*( -COEFB+SQRT ( CQEFB**2-4.*C0EFC I ) /2» 

T4PPOT4P1 
GO TO 75 

90 T2Mi»Ul.-Y2)*T21*Y524T5l«-Y62*ll.-E)*T51)/<l.-Y62*E) 

T3l»EE*T5i«-< l.-EE)*T2Mi 
T61«E*T2MU<i.-E)*T5i 
IF1M.EQ.0I GO TO 51 
M*0 

TT 1PM 1*0.0 

IF1KASE.EQ.1.0R.KASE.EQ.3.0R.KASE.EQ.5) T1PM1«I 1.-YT0T)*T1P1*Y51P4 
IT5UY61P*T61 
IF1T21.EQ.1.0) GO TO 94 
IF1T1P1.EQ.1.0) GO TO 900 
SOI T4P1-T41 1 / I T51-T4PP1 > 

IFISTT.LT*-.OOl) CALL ZERO 12.0, ST, SC, IN,. 0002) 

IF1IN.EQ.4) GO TO 400 

IF! STT.LT.-.OOl.AND. IN.NE.3) GO TO 5i 

94 ETANUNsI l»— YTOT)Pf 1*-T1PI)4| 1«-Y2)*IT1PP1-T2i)-T4P14T41— T51+T4PPI 

95 ETACY»ETANUM/((U-YT0T)*ll.-T61H'(t«-Y2)*tTlPPl-TlPMl)) 
IF(ETANUM.NE.O.O) GO TO 96 

WPARAM*0,0 
GO TO 97 
900 SOO.O 
GO TO 94 

96 WPARAM«« 9487/ETANUH 

97 PllPaPIPlX**<-l./EXP> 

PIPP2*P24PPX**(-l./EXP) 

P4P4-P4P4X**(1./EXP) 

P54PP»P54PPX**(l./EXP) 

IFIJ.NE.O) GO TO 100 



190 IF(J.NE.O) WRITE (6,1060) 

TS1-TS/T1 

IF!SC.LT.O.O.OR.ETACY.LE.O.O.OR.T41.LE.TS1.0R.T4PP1.LE.TS1)GO TO 5 
CAl I RAOTRAI T 1 »TS* T31 » T41 » TAPI * TAPPl , EPS, HR, WPARAM.OLM IN* AINOP, 

IAP ! »,ATOPI 

200 WRITE (6* 1065) T21»E TACY* WPARAH* TIPI *T31 ,T4Pl»T5i *T61*P1 IP*PIPP2*P4P 
l A. P5APP* SC, AROP. AINOP. ATOP 
GO TO 180 

5 WRITEI6* 1070IT21, ETACY, WPARAH, T1P1.T31 , TAPI, T51 ,T6l ,P1 1P.PIPP2, PAP 
HA, P54PP.SC 

180 IF(J.NE.O) GO TO 230 
98 T21-T2UT21DEL 

IFIXT1P1.LE.Y) T21BEG-T21 

IF! ETACY.LT.O.O.AND.T2l~T21BEG.LT. I .1+T21DEL) T21BEG-T21 
I FI ETACY.GE. 0.0) T21END-T21-T21DEL 
IFIT21.LT. T21MAX) GO TO 49 
J«l 

T21-T21BEG 
(NO-1 
GO TO A9 
100 CONTINUE 

CALL HAXIH I T21, T21END.ETACY.I ND..005) 

IFI IN0.NE.6) GO TO A9 
GO TO 190 
230 J-0 

T21-T21NIN 
T4l*T4HT41DEL 
IFI TAI.LT .TA1HAX) GO TO 55 
GO TO (1,2,3), KREAO 
400 (FU.NE.O) GO TO 230 
WRITE (6,2060) T21 
GO TO 180 
999 FORMAT (10(1) 

1000 FORMAT (5F10.A) 

1010 FORMAT! 1HI, A1X,A7HBRAYT0N CYCLE - TWO SPOOL - VARIABLE WORK SPLIT) 

1021 FORMAT! 1HO,AOX,49H1NTERCOOL ONLY - ALTERNATOR ON LOW PRESSURE SMAF 
IT) 

1022 FORMAT! 1 HO, 3 7X,56H INTERCOOL AND REHEAT - ALTERNATOR ON HIGH PRESSU 
IRE SHAFT) 

1023 FORMAT! 1H0,37X,55HINTERC00L ANO REHEAT - ALTERNATOR ON LON PRESSUR 
IE SHAFT) 

1024 FORMAT (1HG,36X,57HN0 INTERCOOL OR REHEAT - ALTERNATOR ON LOW PRESS 
1URE SHAFT) 

1025 FORMAT (1H0,A2X,A6HREHE AT ONLY - ALTERNATOR ON LON PRESSURE SHAFT! 
1030 FORMAT! 1H0,2X,7HETAI Tl ) ,3X, 7HETAIT2) ,3X,7HETA( Cl ) ,3X, THETA! C2) , 5X, 

IAHR! L 1 , 8X , IHE, 7X, 5HGAMMA ,6X,2HC 1 , 8X ,2HCA ,5 X,8H SPLI T IT ) / IOF 10 .4 ) 
1060 FORMAT 1 1H ) 

1035 FORMAT! lH0t5X,2HTl,8X,2HTS,9X,3HEP$,7X,2HHR,l7X,3HYlP,SX,6HVlPFR6t 
16X, 2HY 2, 7X , 5HY 2FR6/AF 10.2* 1OX.4F10. A) 

1045 F0RHATUHK,8H TA/T1-,F6.A,5X,8HT4PP/TI«,F6.A,3X,6HTIPP/TI« ,F6.A, 
l //131H 

1 T2/TI ETACV WCT/P T1P/T1 T3/T1 T6P/TI T5/TI T6/T1 

2 PI/PIP P1PP/P2 PAP/PA P5/PAPP SPLITIC) ARAO/P AINT/P ATOT/P I 
1065 FORMAT! IX, 2F8.A,F9.3,10F8.A,3F8.2) 

1070 FORMAT! IX, 2F8.A,F9.3,10F8. A, 3X,12HN0T COMPUTEO) 

2060 FORMAT! IX, F8. A, 10X.17HN0 VALID SOLUTION) 

ENO 


U ' 
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APPENDIX B 

SUBROUTINE RADTRA 


Subroutine RADTRA computes the specific prime radiator areas required for both 
primary heat rejection and intercooling. This subroutine is entered after cycle tempera- 
tures, cycle efficiency, and weight flow parameter are computed for each combination of 
the cycle temperature variables T41 and T21. It is bypassed during the cycle efficiency 
maximization iteration until the optimum value of T21 is determined. The equations used 
for calculating specific prime radiator area are those in references 1 and 2. 


PROGRAM VARIABLES 

The variable names common to both subroutine RADTRA and main program ETACY2 
are listed in appendix A and omitted here. The variables for RADTRA are the following: 


B 

function of TW3 

CC 

function of EPS and HR 

DELTW 

iteration correction to TW3 

BF 

derivative of B with respect to'TW3 

G 

radiator-area equation term 

LOOP 

iteration counter 

P 

radiator-area equation term 

TW3 

wall temperature at radiator inlet 

TW9 

wall temperature at intercooler inlet 

TW10 

wall temperature at intercooler exit 

TW40 

wall temperature at radiator exit 

T30 

radiator-inlet temperature 

T40 

radiator-exit temperature 

T90 

intercooler-inlet temperature 

T91 

T4P1 in appendix A 

T100 

intercooler-exit temperature 

T101 

T4PP1 in appendix A 
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W radiator- area equation term 

WCPDP function of WCPTOP 

WCPTOP WPARAM in appendix A 


PROGRAM LISTING 


SIBFTC RADTRA 

SUBROUTINE RAOTRAI Tl.TS, T3I ,T4l ,T91 , T101 .EPS, HR, WCPTOP, DLMIN, 
1A I NOP , AROP, ATOP ) 

T30 * T31*Ti 

T40 ■ T4i*TI 

T90 » T91*Tl 

T 100 » T101*T1 

WCP DP-60. 0*WCPT0P/TI 

TM3 - T30 

LOOP * 0 

CC - • 173E-8*EPS/HR 

90 B » TM3-T30«-CC*ITM3**4- TS**4I 
OF * I.0«-4.0*CC*TM3**3 

DELTH » B/OF 

IF ( ABS( 0ELTM/TM3I-DLMINI 91,91,92 
92 TH3 - TM3— OEL T W 
LOOP * LOOP ♦ I 
IF I LOOP-751 90,72,72 

91 TM40 * T40 
LOOP * 0 

95 B » TW40-T40 «-CC*l TM40**4- TS**4I 
OF « 1.0 *4.0*CC*TM40**3 
OELTM - B/OF 

IF I ABSf 0ELTW/TM40I-0LHIN) 96,96,97 
97 TM40 * TM40 - OELTM 
LOOP » LOOP * I 
IF I LOOP-751 95,73,73 
96 I FITS I 9B, 90,99 

98 AROP - NCPDP*f249.0*AL0G<TM3/TM40)/HR420.0*tl.0/TM40**3-I.0 
l/TN3**3)/ I . 1 73E- 8*EF$ ) I 
GO TO 76 

99 P « 60.0*AL0GI ITM3**4- TS**4> /C TM40**4- TS**4))/HR 

M - ALOGI ITM3- TSJMTW40* TSI/IITM40- TSI 
1 HTNimm 

G * 2.0*1 ATANITM3/TS I-4TANITM40/TS I) 

AROP » MCPOP* IP 415.0* I M-G 1/ 1 . 173E-8*EP$* TS**3t) 

76 CONTINUE 
GO TO 50 

72 AROP » 0.0 

GO TO SO 

73 AROP « 0.0 
50 TM9 » TOO 

LOOP • 0 
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54 8 = TH9-T90+CC*(TW9**4- TS**4) 

OF « 1.0*4. 0*CC*TW9**3 
OELTH *• S/OF 

IF ( AB$( DELTH/TW9 I-OLMINI 55,55,56 

56 TH9 ■ TH9-0ELTW 
LOOP » LOOP ♦ I 

IF I LOOP-751 54,70,70 
55 TH 10 * T 100 
LOOP • 0 

57 0 * THIQ-T100*CC*(TH10**4- TS**4) 

of = i,o ♦4.o*cc*rmo4*3 

OELTH » 3/OF 

IF ( ABSI DELTW/THIOI-DLMINI 59,59,60 
60 THIO » THIO -OELTH 
LOOP * LOOP ♦ I 
IF (LOOP-751 57,71,71 
59 IF( TS 162, 62, 63 

62 A (NOP * WCPDP*(240.0*AL0G( TH9/TW10) /HR*20. 0*11 
l/IH9**3l/( .1 73E-8*EPS 1 1 

GO TO 51 

63 P « 60.0*ALOG( (TH9**4- TS**4) /( TH10**4- 

H * AL0GKTH9- TS)*(TH10* TSI/((TH10- 

1 *<TH9*TS))I 

G * 2»0*( ATANl TH9/TS I-ATANCTHIO/TS it 

AINOP * WCPQP*(P*15.»04( H-G) / ( • 173E-8*EP$4 
GO TO 51 

70 AINOP » 0.0 

GO TO 51 

71 AINOP * 0.0 

51 ATOP * AINOP ♦ AftOP 
RETURN 
END 


• 0/ T Hl0**3-1.-*0 

TS**4I I /HR 
TS) 

TS*93I) 



APPENDIX C 

SUBROUTINE MAXIM 


Subroutine MAXIM provides the logic for determining the maximum cycle efficiency 
and associated turbine temperature ratio (T21) within the range specified by T21BEG and 
T21END. For each cycle temperature ratio, this subroutine is entered after the para- 
metric results for all turbine temperature ratios have been determined. The subroutine 
logic causes the maximum function value to be approached from both sides until there is 
a convergence at the peak. 


PROGRAM VARIABLES 

The variables (or MAXIM are as follows: 


JUMP 

logic indicator for branching 

SPEED 

retained values of WA 

TOLER 

convergence tolerance 

WA 

T21 in appendix A 

WAMAX 

T21END in appendix A 

WEIGHT 

retained values of WTFL 

WTFL 

ETACY in appendix A 


PROGRAM LISTING 


SI8FTC NAXIN 

SUBROUTINE NAXIN fWA v WANAX«WTFL , INO, TOLER) 

01 MENS ION SPEED! 3) •WEIGHT! 3) 

C INO * I v FIRST POINT 
C INO • 2, SECOND POINT 

C INO « 3. THIRD POINT (EUNINATE I POINT IF NECESSARY) 
C INO « 4» FOURTH POINT f ELININATE 4TH POINT) 

C INO - 5* WANAX IS LESS THAN HA 

GO TO <l40 v l50 v 2t0 9 270 ( 370)«IND 
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non 


140 JUMP * l 

IF ( WAMAX.LT. WAJ GO TO 145 
SPEEOm * WA 
WEIGHT ( 1 ) * WTFL 
WA * WANAX 
INO * 2 
RETURN 

145 SPEED! 3) * WA 
WEIGHT(3I * WTFL 
WA * WANAX 
INO * 5 
RETURN 

150 SPEEOf 3) » WA 
WEIGHT! 3) * WTFL 
160 WA a ! SPEED! 1 ) ♦SPEED! 3) )/2« 

IF! I SPEEO! 31-SPEED! 1 I ) .LT.TOLER) GO TO 400 
INO * 3 
RETURN 

210 SPEEO! 21 > WA 
WEIGHT! 21 « WTFL 

I FIWTFL.LE .WEIGHT! ll.OR.WTFL.LE.WEI GHTI3) ) GO TO 268 
INO > 4 

CHOOSE PROPER INTERVAL FOR NEXT POINT !4TH POINT HAS BEEN ELIHINATEOI 

245 GO TO !246«247)t JUMP 

246 JUMP * 2 
GO TO 250 

247 JUMP a l 
GO TO 260 

250 WA * (SPEEO! 1) ♦SPEED! 2) 1/2.0 
RETURN 

260 WA * | SPEED! 31 ♦SPEEO! 211 /2.0 
RETURN 

268 IF! WEIGHT! 3) .GT. WEIGHT! 1)1 GO TO 269 
WEIGHT! 3) « WTFL 
SPEEO! 31 « WA 
GO TO 160 

269 WEIGHTU) a WTFL 
SPEED! 1) a WA 
GO TO 160 

270 IF! ( SPEEO! 31-SPEED! 1 I I.LT.TOLER) GO TO 400 
280 IF (WTFL ■'HEIGHT! 2) I 320.350,290 

C 

C NEW POINT BECOMES NIDPOINT v MIDPOINT BECOMES END POINT 
C 

290 IF (WA-SPEE0(2I) 310*300*300 
300 SPEEO! II « SPEEO! 21 
SPEEO! 21 a WA 
WEIGHT! II a WEIGHT (2) 

WEIGHT! 2) a WTFL 
GO TO 295 

310 SPEED! 31 a SPEED! 21 
SPEED! 21 a WA 
WEIGHTI3I • WEIGHT! 21 
WEIGHTI2I « WTFL 
GO TO 265 



non 


NEW POINT BECOMES END POINT 

320 IF ! WA-SPEED! 2 ) I 340,330*330 

330 WEIGHTI3I = WTFL 
SPEED! 31 * WA 
GO TO 245 

340 WEIGHT! 1 1 = WTFL 
SPEED! 1) « WA 
GO TO 245 

350 IF! WA«GTcSPEED(2) ) GO TO 360 
SPEED! 3) * SPEED! 21 
WEIGHT13I * WE IGHT! 2) 

GO TO 210 

360 SPEED! 1) * SPEED! 2) 

WE IGHT! 1) = WE IGHT! 21 
GO TO 210 

370 SPEEDil) - WA 

WEIGHT! 1 ) * WTFL 
GO TO 160 

400 INO * 6 
RETURN 
ENO 
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APPENDIX D 
SUBROUTINE ZERO 


Subroutine ZERO provides the logic for determining the turbine work split (ST) that 
yields a two- shaft one- compressor arrangement (SC = 0). When the input value of ST is 
negative, subroutine ZERO is entered after each calculation of SC in order to test the 
value of SC and select a new estimation for ST. A linear interpolation or extrapolation is 
used to determine the estimated value of ST. 


PROGRAM VARIABLES 

The variables for ZERO are as follows: 

TOL convergence tolerance 

X ST in appendix A 

XX value of ST for second iteration 
XY reciprocal of slope of estimation line 

XI X value of first point defining estimation line 

X2 X value of second point defining estimation line 
Y SC in appendix A 

Y1 Y value of first point defining estimation line 
Y2 Y value of second point defining estimation line 

PROGRAM LISTING 


SIBFTC ZERO 
C 

C FOR NONOTONIC DECREASING Y*FCN(X) 
C 

SUBROUTINE ZERO IXX,X. Y, IN, TOL) 
IF(ABSfYt.LE.TOL) GO TO 100 
GO TO tlOt 20)* IN 


23 



10 Xl*X 
Y l*Y 
x=xx 
IN*2 
RETURN 
20 X2*X 
Y2=Y 

XY»(X2-X1»/(Y2-Y1> 

IF< XY.GF,0.0) GO TO 200 
X*XY*f-Yll+Xl 

IF( ABS(Y1)-A8S(Y2) ) 50,40*40 
40 X1*X2 
Y1*Y2 
50 RETURN 
100 IN* 3 

Y*A8S( V) 

RETURN 
200 IN*4 
RETURN 
ENO 
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